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ABSTRACT 

Three-dimensional (3D) spectroscopy techniques are becoming more and more 
popular, producing an increasing number of large data cubes. The challenge of 
extracting information from these cubes requires the development of new tech- 
niques for data processing and analysis. We apply the recently developed tech- 
nique of Principal Component Analysis (PCA) Tomography to a data cube from 
the center of the elliptical galaxy NGC 7097 and show that this technique is effec- 
tive in decomposing the data into physically interpretable information. We find 
that the first five principal components of our data are associated with distinct 
physical characteristics. In particular, we detect a LINER with a weak broad 
component in the Balmer lines. Two images of the LINER are present in our 
data, one seen through a disk of gas and dust, and the other after scattering 
by free electrons and/or dust particles in the ionization cone. Furthermore, we 
extract the spectrum of the LINER, decontaminated from stellar and extended 
nebular emission, using only the technique of PCA Tomography. We anticipate 
that the scattered image has polarized light, due to its scattered nature. 

Subject headings: Galaxies: active - galaxies: elliptical and lenticular, cD - galax- 
ies: individual (NGC 7097) galaxies: kinematics and dynamics - galaxies: nuclei 
- techniques: spectroscopic 



1. INTRODUCTION 

Active Galactic Nuclei (AGN) form a heterogeneous group of objects including quasars, 
Seyfert galaxies, LINERS and others. The dominant view today is that these objects are 
associated with the capture of gas by a supermassive black hole located at the center of 
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those galaxies. In the optical, their spectra are characterized by broad and intense emission 
lines. Quasars are among the more luminous types of AGN and are typically found at higher 
redshifts. At lower luminosities one finds Seyfert galaxies, whose spectra are characterized 
by intens e lines from hig h ionization species and LINERs (Low Ionization Nuclear Emitting 
Regions - lHeckmanlll980l ) . with enhanced lines from low ionization species. In both cases, the 
host galaxy is usually detected. Seyfert galaxies are typica lly spiral galaxies, while LINERs 
and quasars are more frequently seen in elliptical galaxies (jHol 120081 ). 



Seyfert galaxies (and LINERs as well) can be divided into two types: those that have 
permitted lines wider than the forbidden lines (type 1) and those that have both types o f 



lines with similar widths (type 2). The unified model (lOsterbrockl Il978l ; lAntonucci 



19931 ) 



suggests that the broad line emission region (BLR) is surrounded by a torus of gas and 
dust; so depending on the orientation of this torus relative to our line of sight, one can or 
not observe the BLR, thereby explaining the existence of the two types of Seyfert galaxies. 
This model is convincingl y verified by the spectropolarimetric observations of NGC 1068 
(lAntonucci k Miller! E9851 ). 



In this article, we analyze a data cube for the central part of the elliptical galaxy NGC 
7097, obtained with the GMOS-IFU spectrograph on the Gemini South Telescope, with 
the goal of detecting and characterizing low-luminosity AGNs in early type galaxies. Data 
cubes are large sets of measurements with two spatial (spaxels) and one spectral (pixels) 
dimension, composed of several tens of millions of measurements. Extracting information 
from such cubes is a complex task and the traditional methods of spectral analysis are 
inadequate, so new methods are much needed and welcome. The number of spectrographs 
capable of producing such large data cubes is increasing rapidly. For instance, the JWST 
(James Webb Space Telescope) will have 5 IFUs and all projected ELTs (Extremely Large 
Telescopes) plan to provide at least one IFU in their first generation instruments. In this 
paper, we an alyze a data cube o f NGC 7097 with a recently developed methodology (PCA 
Tomography; ISteiner et al.l 120091 ) and show its power to extract useful information. 



NGC 7097 is an E5 ell iptical galaxy, acc ording to RC 3 (Ide Vaucouleurs et al.l Il99ll ) 
at a distanc e of 32.4 Mpc (ITonry et al.l 120011 ). This galaxy shows strong emission of [O 
II] A3727A Jcaldwelj ll984j) and [N II] A6584A which has a widt h of FWHM = 423 km/ s 
( Phillips et al. 19861 ); the strong [N II] emission was confirmed by Pizzella et al. ( 1997 ). The 
presence of these strong low ionization emission lines, not produced in H II r egions, indicates 
that the gas requires a non-thermal ionization source. ICaldwell et al.l (Il986[ ) detected a disk 
of gas, with a radius of 15" , aligned with the axis of the galaxy (the position angle of the 
galaxy is PA = 18°), counter-rotating with respect to the stellar component. As the extent 
of the galaxy is 1.9 arcmin on the sky, the gaseous disk has nearly \ of the size of the galaxy. 
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However, iBuson et al.l (119931 ). using imaging in Ha+[N II] and the adjacent continuum, 
found that the extended emission is rotated by —30° from the stellar isophotes, that is, the 
emissio n disk has a PA = — 12°. These authors estimate the inclination of the gaseous disk 
as 57°. IZeilinger et al.l ( 119961 ) also de tected a disk s t ructu re with an inclination i = 62° with 
respect to the plane of the sky while iPizzella et al.l (Il997| ) found 67°. 



2. OBSERVATIONS, REDUCTION AND DATA TREATMENT 

The observations of NGC 7097 were made on July 31, 2007 at the Gemini South Tele- 
scope 0. We used the "Integral Field Unit" - IFU of the Gemini Multi Object Spectrograph 



GMOS (IHook et al. 2004 



Allington- Smith et al.ll2002l ) in the single slit mode. This mode 



allows the object and the sky to be observed simultaneously by 500 and 250 micro-lenses, 
respectively, located in the focal plane of the telescope. The micro-lenses divide the image of 
the object into slices of 0.2 arcseconds and is coupled to an array of optical fibers arranged 
linearly on the nominal location of the slit spectrograph (pseudo-gap). The observation of 
the object and sky are separated by one arcminute. The final product is a data cube with 
two spatial dimensions spanning 3.5 x 5 arcsec and one spectral dimension. This allows 
the construction of images in a specific wavelength range, or the extraction of the spectrum 
in a given region of space. We use the grating B600-G5323, with a central wavelength of 
observation in 5650A. The spectra cover a range of 4228 - 7120A and have a resolution of 
1.8 A, as measured from the skyline [O I] 5577A. 

The usual lamp and sky flat-fields as well as bias exposures were taken for the data cube 
corrections. Exposures of a Cu-Ar lamp were taken for wavelength calibration and a data 
cube of the star LTT 9239 was obtained for the flux calibration. The seeing was measured 
to be 1.0". The resulting data were reduced using the standard Gemini IRAF0 package and 
a spatial sampling of 0.05 arcsec per pixel, while the CCD spaxel has 0.1" x 0.1". Bias and 
background subtractions were made, as well as the corrections for the fiber response, removal 



1 Based on observations obtained at the Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of 
the Gemini partnership: the National Science Foundation (United States), the Science and Technology 
Facilities Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Australia), Ministerio da Ciencia e Tecnologia (Brazil) and Ministerio de 
Ciencia, Tecnologia e Innovation Productiva (Argentina) 

2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Asso- 
ciation of Universities for Research in Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 
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of cosmic rays and subtraction of the sky from both the observations of the galaxy and also 
from the standard star. 

T he spatial high frequency noise in the data cube was suppressed using a Butterworth 
filter ( Gonzalez fc Woods! 120021 ) applied to the Fourier transform of each image in the data 
cube. The parameters of the filter H(u,v) were n = 6 and a = b = 0. 15 Fny, i.e. H(a,b) = 
0.5, where Fny is the Nyquist frequency of the spatial component of the cube. After removing 
the high frequency noise, the data cube of NGC 7097 still has a signature of low-frequency 
instrumental noise, which is present in almost all the observations made with the GMOS- 
IFU. To reduce such instrumental "fingerprints" we have developed a technique that uses 
PCA Tomography in the wavelet space of the data cube's spatial dimensions. We have no 
space here to describe this procedure, nor is it the scope of this letter. This technique is 
laborious and is described in detail in Steiner et al. (2011, in preparation). This technique 
was applied to successfully remove the low frequency instrumental noise to the data cube 
discussed in this Letter. 

The GMOS-IFU data were obtained without an atmospheric dispersion corrector (ADC). 
Thus, the data cube suffers from the effect of atmospheric dispersion, which increases with 
the air mas s of the observati on. To correct for the atmosp heric dispersion, we used the 
equations of iFilippenkol ( 119821 ) and lBonsch fc Potulskil ( 119981 ). Finally, the cube was decon- 
volved in the spatial dimension using the Richardson- Lucy method with 6 iterations and 
a Gaussian PSF with FWHM = 1.0". After deconvolution, the estimated final PSF has a 
FWHM = 0.7", as indicated by our experience from processing many other stellar profiles 
of other observing programs. 



3. DATA ANALYSIS AND RESULTS 

The analysis of large data cubes may become complex and overwhelming, as it may 
involve tens of millions of pixels. More concerning is that, given this complexity, only some 
restricted subset of the data ends up being analyzed; the rest is at the risk of being largely 
ignored. To overcome the complexity of analyzing such large data cubes, we employ a method 
of data cube analysis based on Principal Component Analysis (PCA). This method extracts 
the significant information content associated with the data through an effective dimensional 
reduction, facilitating its interpretation. PCA compresses the data, originally expressed as a 
large set of correlated variables, into a small but optimal set of uncorrelated variables, ordered 
by their eigenvalues. An important aspect of PCA is that the eigenvectors are mutually 
orthogonal; if there are physically uncorrelated phenomena expressed in the data cube, they 
possibly will align with different eigenvectors in decreasing order of importance. Furthermore, 
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PCA is a nonparametric statistic, therefore there are no parameters or coefficients to adjust 
that somehow depend on the users' experience and skills, or on physical and geometrical 
parameters of a proposed model. PCA, therefore, provides a unique and objective solution. 

PCA has been used many times in t he astronomical lit e ratur e, an d a more e x tende d 
presentation of this technique is given in iMurtagh fc Heckl (119871 ) and iFukunagal ( 11990I ). 
Most of the applications of PCA in astronomy are related to finding eigenvectors for a 
population of objects. In the present case, we want to apply the technique to a single 
data cube in which the objects are spaxels (spatial pixels). The wavelength pixels are 
their properties. PCA produces eigenvectors (the uncorrelated variables), which we refer 
to, also, as eigenspectra, and tomograms, which are images of the data projected onto each 
eigenvector. In traditional tomographic techniques, one obtains images that represent 'slices' 
in three-dimensional space (the human body, for example), or in velocity space (Doppler 
Tomography). In PCA Tomography, one obtains images that represent 'slices' of the data 
(tomograms) in the eigenvector space. Each tomogram is associated with an eigenspectrum. 
The simultaneous analysis of both brings a new perspective to the in terpretation of them. 
For a full pr esentation of PCA Tomography, see ISteiner et al.l (120091 ); a shorter version is 
presented in ISteiner et al.l ( 120101 ) . 



We will illustrate the method of PCA Tomography by applying it to our data of the 
galaxy NGC 7097 and interpreting the first 5 principal components in terms of physical 
properties. The first eigenvector, which explains 99.53% of the data cube variance is shown 
in Figured] and is very similar to the integrated spectrum of the field of view (see Figure 
The added contribution of eigenvectors 2-5 explain 0.46% of the variance in the data cube. 
This shows the great redundancy of this type of data and has relevance in discriminating 
non-redundant information. By removing the effects of the strongest correlations one can 
look to the less significant ones. The tomograms and eigenvectors and 2-5 are shown in 
Figure [31 

Eigenvector 2, explaining 0.38% of the variance, reveals correlations between features 
associated with transitions of H/3, [O III] AA4959A, 5007A, [O I] A6300A, Ha [N Hi AA6548A , 
6583A and [S II] AA6714A, 6732A, with relative intensity typical of LINERs JHeckmanlll98ol ). 
These characteristics are also correlated with a red continuum and interstellar absorption of 
Na I A5891A. The tomo gram shows the location of the object. 

In eigenvector 3, explaining 0.05%, we observe an anti-correlation between the red and 
blue wings of the features associated to the emission lines (also see Figure H]). The tomogram 
reveals that the blue (approaching) material is in the south while the red (receding) material 
is in the north. The position angle (P A) is —10°. This is a signature of rotating gas, which 



is also seen in the galaxies NGC 4736 (ISteiner et al.ll2009l ) and M 81 (jSchnorr Miiller et al. 
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Fig. 1. — Eigenvector (Ev 1) and tomogram of principal component 1 of the data cube of 
NGC 7097. Black means stronger and white, weaker correlation 
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Fig. 2. — Spectrum of NGC 7097. The top curve shows the spectrum summed over all 
spaxels. This spectrum is dominated by the stellar emission in the field of view, but we can 
see the lines of Ha, [N II] and [S II]. The bottom curve shows the extracted AGN spectrum, 
using the technique of feature suppression. In this spectrum we can clearly see the features 
of LINER emission and broad, although weak, Balmer lines. The two spectra are shown in 
relative intensity units, but on the same scale. 



20111). 
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In the fourth eigenvector, with 0.02%, of the variance, there is again a correlation be- 
tween emission line features typical of a LINER, but this time correlated with a blue contin- 
uum and the absence of interstellar absorption of Na I. The radial velocities in eigenvectors 
2 and 4, estimated from the [O I] A6300A transition, differ by less than 10 km/s, well below 
the spectral resolution, which is 90 km/s. 

The fifth eigenvector, with 0.01% of the variance presents an eigenspectrum in which 
there is an anti-correlation between the low-velocity line features and their blue and red 
wings (see Figures [3] and H]) . The fifth tomogram reveals a low- velocity extended correlation 
that is located along the axis of rotation of the gas (perpendicular to the disk). The region 
associated with the line wings, for the high rotation velocity, corresponds to the disk, as seen 
in the third eigenvector. We interpret the low velocity region as the ionization bi-cone. 

How far can we look to the eigenvector befo re running into noise? We have made the 
"scree test" (see Figure 1 in ISteiner et al.l 120091 for an explanation) and found that noise 
dominates all eigenvectors above number 7. Eigenvalue 8 is 0.0017%, 10 times smaller than 
eigenvalue 5. 

The LINER appears twice, in components 2 and 4. It is possible to reconstruct a data 
cube removing eige nvector 1 using the technique of "feature suppression" (see section 6 of 



Steiner et al 



20091 ). Performing the suppression allows us to extract the spectrum of the 



central region, shown in Figure [2J This illustrates, using only PCA, how one can isolate the 
spectrum of a low luminosity AGN decontaminated from the stellar and nebular emission. 
In Figure [2] we see that the Balmer lines have a weak broad component. This object is, 
therefore, of type 1 (sometimes also classified as type 1.8, given the weakness of the broad 
Ha, o r Lb). The relative intensities of the narrow lines are consistent with LINER (IHeckman 
198(j): we measured [O III] A5009A / Hff ~ 1.5, [N II] A6584A / Ha ~ 1.2 (compare with 
Pizzella et aD Il997l and IPhillips et al.l Il986h and [O I] A6300A / Ha ~ 0.34. For typical 
LINERs, [O IHl A5009A / Hff < 3, [N II] A6584A / Ha > 0.5 and [O I] A6300A / Ha > 0.05 
( IQsterbrock fc Ferlandl 120061 ) . From the original spectrum, we estimate the reddened Ha 
luminosity as 5.0 ± 1.0 x 10 37 erg s _1 . We will not attempt to estimate the reddening as this 
will be done in more detail in a forthcoming paper. 



4. DISCUSSION AND CONCLUSIONS 

We have seen that the eigenvectors and tomograms 2, 3, 4 and 5 can be interpreted 
as physical phenomena with distinct spectral properties. We summarize tomograms 2-5 in 
Figure [5J On the left we see each one with a specific color (as indicated) and on the right we 
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show a sketch that illustrates our interpretation. It seems clear that the LINER is seen twice: 
once directly through the disk of gas and dust (Ev 2), manifested by the heavily extincted 
continuum spectral slope and by interstellar neutral gas absorption. This same AGN is also 
seen (Ev 4) scattered by free electrons existing in the ionization cone; dust particles could 
also be present and produce additional scattering, especially in the blue. As we do not have 
reliable information on the absolute spectral behavior of the reflected continuum, we cannot 
say whether the scattering is independent of wavelength (as expected by electron scattering) 
or blue (if the scattering is produced by dust). This scattered image is less subjected to 
dust reddening or interstellar absorption by the disk, at least not as much as the object seen 
directly. 

The gaseous disk, identified in ei genvector 3, has a pos ition angle of PA = —1 0°; it is 



co-aligned with the disk described by ICaldwell et al.l (119861 ) and iBuson et al.l ( 119931 ) which 
is ~ 15 times larger. 

We identify eigenvector 5 as representing the bi-cone in anti-correlation with the disk. 
The profile has a correlation in the low velocity core anticorrelated with the wings. Spatially, 
the wings are associated with the region corresponding to the disk (EV3). The central 
core in the eigenspectrum (see Figures |3] and H]) , seen in anti-correlation with the wings, is 
correlated with a structure in the respective tomogram that is perpendicular to the gaseous 
disk-looking, just like an ionization bi-cone. This disk-cone anti-correlation is expressed 
because the disk is quite highly inclined, presenting high velocities while the cones are less 
inclined with low intrinsic velocities. In agreement with this interpretation, an image of 
the [O III]A5007A shows a weak extended emission along the near-side cone. One problem 
with this interpretation is that the far-side cone is slightly brighter than the near-side. We 
have to assume that this is an intrinsic asymmetry. Although this is not necessarily the only 
possible interpretation for the eigenvector 5, it seems to be self-consistent and is the only one 
we could find. The existence of this ionization bi-cone has not been reported previously for 
this galaxy. In drawing the d i sk sh own in Figure |5l we have assumed the inclination angle 



determined by iPizzella et al.l ( 119971 ). 67°, and postulated a cone half-opening angle of 57°. 
An opening angle much different from this would not produce two such close AGN images 
as observed. 

The separation between the two images of the AGN is about 0.2" while the seeing was 
~ 1" (Figure [5]). This is possible because the PCA is a filtering technique. Each tomogram 
produces an image of correlations and, therefore, can be interpreted representing the data 
in distinct "filters" . This does not produce super-resolution but can image phenomena with 
angular separations smaller than the resolution. The basic reason for this is that correlations 
are present below the seeing limit and not destroyed by the random variations in the PSF. 



-9 - 



We predict that the blue dot in Figure [5] should have polarized light because of its scattered 
nature. 

In conclusion, we presented a new method of locating and extracting scattered light 
from AGN; this does not use the well known method of spectropolarimetry but uses PCA 
analysis combined with 3D spectroscopy. 

We would like to thank Alex Carciofi and J. E. Bjorkman for critically reading the 
manuscript. We also thank Felipe Andrade Oliveira for helping in the ilustration. T.V.R. 
and R.B.M. thank FAPESP for support under grants 2008/06988-0 (T.V.R) and 2008/11087- 
1 (R.B.M.) . We would like to thank an anonymous referee for his valuable suggestions that 
improved the quality of this Letter. 
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Fig. 3. — Eigenvectors (Ev2 to Ev5) and tomograms of principal components 2-5 of the data 
cube of NGC 7097. In tomograms 3 and 5, "black" means correlation and "white" means 
anticorrelation. In tomograms 2 and 4, we display only the tip of the positive correlations. 



-12- 





4980 4990 5000 5010 5020 5030 
Wavelength (A) 



Fig. 4. — The expanded eigenvectors 3 (Ev 3 - top) and 5 (Ev 5 - bottom) for the wavelength 
region near the [O III]5007A line. These characteristics are interpreted as the anti-correlation 
produced by a rotating disk (in eigenvector 3) and that between the disk (wings of the lines) 
and the bi-cones (core of the lines) (in eigenvector 5). 





Arcsec 

Fig. 5. — Tomograms 2, 3, 4 and 5. The left image shows the combination of tomograms 2 
(magenta), 3 (red=positive; blue=negative), 4 (cyan) and 5 (green). The extended emission 
in blue and red corresponds to the rotating disk of gas around the AGN. The small dots 
correspond to the emission of the AGN seen directly (magenta) and reflected in the ionization 
cone (cyan). The emission in green represents the bi-cone ionization. Right: Model of the 
central region of the galaxy NGC 7097. 



